Transport properties of the metallic state of overdoped cuprate superconductors from an 

anisotropic marginal Fermi liquid model 
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We consider the implications of a phenomenological model self-energy for the charge transport properties 
of the metallic phase of the overdoped cuprate superconductors. The self-energy is the sum of two terms with 
characteristic dependencies on temperature, frequency, location on the Fermi surface, and doping. The first 
term is isotropic over the Fermi surface, independent of doping, and has the frequency and temperature de- 
pendence characteristic of a Fermi liquid. The second term is anisotropic over the Fermi surface (vanishing 
at the same points as the superconducting energy gap), strongly varies with doping (scaling roughly with T c , 
the superconducting transition temperature), and has the frequency and temperature dependence characteristic 
of a marginal Fermi liquid. Previously it has been shown this self-energy can describe a range of experimen- 
tal data including angle-dependent magnetoresistance (ADMR) and quasi-particle renormalisations determined 
from specific heat, quantum oscillations, and angle-resolved photo-emission spectroscopy (ARPES). Without 
introducing new parameters and neglecting vertex corrections we show that this model self-energy can give a 
quantitative description of the temperature and doping dependence of a range of reported transport properties of 
T12201 samples. These include the intra-layer resistivity, the frequency dependent optical conductivity, and the 
Hall coefficient. The temperature dependence of the latter is shown to be particularly sensitive to the anisotropy 
of the scattering rate and to the shape of the Fermi surface. In contrast, the temperature dependence of the Hall 
angle is dominated by the Fermi liquid contribution to the self-energy that determines the scattering rate in the 
nodal regions of the Fermi surface. 

PACS numbers: 74.72.-h, 74.72.Gh, 74.62.-c, 75.47.-m 



I. INTRODUCTION 

Much research of strongly correlated electron materials, 
with high-temperature superconducting cuprates being the 
prominent example, is focused on the experimental or theoret- 
ical determination of the relevant electronic self-energy. That 
is because the self-energy can provide insight into the under- 
lying quantum many -body physics. The proper knowledge on 
the self-energy in the metallic phase of high-temperature su- 
perconductors is also believed to be a step towards solving 
the mystery of high-temperature superconductivity since ulti- 
mately superconductivity is an instability in the metallic state. 
A model self-energy capable of a unified description of results 
from many experiments is therefore desirable and provides a 
benchmark for comparison with microscopic theories based 
on lattice effective Hamiltonians such as the t — J and Hub- 
bard models. 

In the last two decades experimental data has been 
used to deduce the self-energy, both directly and indi- 
rectly. Angle-resolved photoemission spectroscopy (ARPES) 
offers information on both the real and imaginary part 
of the self-energy^, specific heat provides information 
through renormalization effects^!, and angle-dependent 
magnetoresistance^— (ADMR) provides information on the 
imaginary part of the self-energy or scattering rate close to 
the Fermi level. Further information about the temperature 
and doping dependence can be obtained from measurements 
of the resistivity 15-19 and Hall effect!^ 2 !. In addition, the op- 
tical conductivity 22-25 provides information on the frequency 
dependence of the self-energy. 

Previous work 26 introduced a particular model self-energy, 



motivated by ADM R 11 ' 12 , that could describe consistently 
and quantitatively ADMR and a number quantities deter- 
mined by the real part of the self-energy, including ARPES 
dispersion 8 , specific heat 910 , and effective masses deduced 
from quantum oscillations 27 , in the entire overdoped regime 
of Tl2Ba2Cu06+5 (T12201). In this paper we extend our anal- 
ysis to description of transport properties, which are largely 
determined by the imaginary part of the self-energy. Prop- 
erties considered include the intra-layer resistivity^ - — , Hall 
effec t 17,18,20 ' 21 , and optical conductivity 22-24 . We show that all 
of these can be quantitatively described with our model self- 
energy without any additional fitting parameters for the entire 
overdoped regime for T12201 . 

The outline of the paper is as follows. In Section II we re- 
view the form of the self-energy and its parametrisation. Sec- 
tion III considers the DC conductivity and the frequency de- 
pendent conductivity. It is shown that at high temperatures 
and frequencies these are sensitive to the cut-off frequency 
which appears in the Fermi liquid term in the self-energy. In 
Section IV we show that the Hall coefficient strongly depends 
on the shape of the Fermi surface and that its non-monotonic 
temperature dependence gives strong support for our model 
self-energy. We also argue that the observed non-monotonic 
temperature dependence of the Hall coefficient cannot be cap- 
tured with some alternative models, e.g., with the isotropic 
marginal Fermi Liquid model 28 . We also argue that the ob- 
served T 2 dependence of the Hall angle arises because it is 
dominated by the isotropic part of self-energy, which also 
equals the smallest scattering rate on the Fermi surface in the 
nodal direction, and that the contribution of the anisotropic 
part is suppressed. Hence, results on the Hall angle give ad- 
ditional support to our model self-energy, in particular the T 2 
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dependence of the isotropic part of self-energy. In Section V 
we briefly review relevant results from microscopic model cal- 
culations. Although, several are qualitatively consistent with 
the model self-energy, they tend to obtain scattering rates that 
are significantly less than observed. Section VI contains some 
conclusions and suggestions for possible future work. 

n. MODEL SELF-ENERGY 

Our model self-energy is motivated by the angle- 
dependent magnetoresistance (ADMR) experiments on over- 
doped T12201— ~—, where two distinct scattering rates were 
uncovered. The first is more Fermi liquid (FL) like and 
is isotropic over the Fermi surface, weakly doping depen- 
dent, and shows T 2 dependence at low T. The second has a 
marginal Fermi liquid 2 ^— frequency and temperature depen- 
dence and is strongly anisotropic over the Fermi surface (the 
same anisotropy as the superconducting gap). Its strength fol- 
lows the doping dependence of T c in the strongly overdoped 
regime and is linear in T down to the lowest T. 

Accordingly, our model self-energy can be written, 

S"(k, W ) = S» + Sl ME (^), (1) 

where <fi denotes the position on the Fermi surface (see Fig. 
01. The imaginary part of the isotropic FL like self-energy is 
given by22 

fi_ a sd±gZ± for " 2 +: 2 / 2 < 1, 

(2) 

Here 1/(2to) accounts for the impurity scattering, and 
Matthiessen's rule is implicitly assumed. The parameter s 
gives the strength of the FL like self-energy part and ui^ is 
the high-w cutoff (see Fig. [TJ. We use units h = ks = 1. SpL 
is quadratic in uj and T at low uj and T. The function F is a 
slowly decreasing function with F(l) = 1, which we simply 
approximate with a constant. 

The anisotropic marginal Fermi liquid (AMFL) part of the 
self-energy (its imaginary part) has the following form, 

where A(<fi) determines its strength and is anisotropic over the 
Fermi surface, x = max(|w|, 7rT) and cl^mfl is the high-w 
cutoff. S^( MFL (0, uj) is linear in to and T for low u or low T 
(see Fig. [TJ. The real part of the self-energy is obtained from 
a Rramers-Kronig relation and is not explicitly given here. 
Explicit low uj behaviour of the real part can be found in Ref. 

m 

Parameters of the model self-energy were already extracted 
in Ref. |26|for overdoped T12201. From ADMR one can esti- 
mate that s/w F l = 9.2 eV" 1 and 

\(4>) = 1.6 cos 2 (20)T c (p)/T c max , (4) 

where T c (p) is the doping (p) dependent transition temper- 
ature and T™ ax is the maximal transition temperature. For 
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Figure 1 . Frequency dependence of the imaginary part of the model 
self-energy. The first part Efl, is Fermi liquid like, with a quadratic oj 
and T dependence, up to the high-o; cutoff ojp L , and with a prefactor 
s/ujft, This part is taken to be doping independent. The second part 
is the anisotropic marginal Fermi liquid part, whose imaginary part 
is constant for oj < nT and proportional to T, while it is linear in 
ui for higher uj up to a high-o; cutoff oj^mfl- I ts strength is given by 
A(0), which is strongly anisotropic over the Fermi surface and also 
strongly doping dependent. 



the doping dependence of T c we use the phenomenological 
relation 3 ? -T c (p)/T™ x = 1 - 82.6(p - 0.16) 2 with T c max = 93 
K at the optimal doping p = 0.16 for T12201. This T c (p) 
relation is for illustrative purposes only since the supercon- 
ductivity actually survives up to p = 0.31 as was found in 
Ref. |34j In addition, Wp L was estimated 26 to be 0.23 eV from 
specific heat measurements in the highly overdoped and non- 
superconducting regime. The cutoff w^mfl on ty weakly in- 
fluences the results 26 because it only enters the real part of the 
self energy via a logarithmic dependence and is here taken to 
be 0.2 eV. For the T12201 samples used in ADMRii^ the im- 
purity scattering rate was estimated to be, 1/(2to) ~ 4 meV. 
These parameter values for the model self-energy give a con- 
sistent description of several experiments, including ADMR, 
specific heat, quantum oscillations, and the quasi-particle dis- 
persion seen in ARPES^ 6 -. 

For Hubbard models, there is an additional constraint on the 
self-energy, and in particular its high-frequency behaviour, via 
the sum rule^S 

J du(-K(k,u)) = 7r£/n_ CT (l-n_ ff ), (5) 

with U being the on site Coulomb interaction strength, and 
n a = (l+p)/2 being the density of electrons with spin a. Our 
model self-energy does not obey this sum rule since S"(fe, uj) 
stays finite for uj — > oo. To fulfill the sum rule our S"(fc, uj) 
should be strongly suppressed at high frequencies. We esti- 
mate this suppression should occur at uj ~ 5 eV (for U = 8i). 
Such a suppression would not influence our results, since they 
are determined by the value of the self energy at much lower 
frequencies. Hence, we do not employ the suppression and the 
self-energy sum rule in this work. In contrast, in the next Sec- 
tion it is shown that the behaviour of the self-energy near our 
cutoff frequencies u>p L and w^ MFL does affect some observed 
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transport properties at high temperatures and frequencies. 

For the bare band dispersion eo(A:) we approximate the 
LDA results in Ref. 37 with the following hopping param- 
eters, ti = 0.438, t 2 = -0.150, t 3 = 0.084, t A = -0.013, 

= —0.020, t§ = 0.029, all expressed in eV (for details see 
supplemental material of Ref. l2q) . To obtain the Fermi sur- 
face volume of the overdoped regime we apply a rigid band 
shift through the chemical potential /i. The main doping de- 
pendence of our results does not come from the band filling, 
which we therefore keep fixed, but rather from the doping (or 
T c ) dependence of the self-energy. Shifting the chemical po- 
tential from values appropriate for highly overdoped to opti- 
mal doping (e.g., from p = 0.3 to 0.15) induces only a small 
change in the Fermi surface (see Fig. |8]for example). 



III. INTRA-LAYER CONDUCTIVITY 

The frequency dependent conductivity is approximated 
with the bubble diagram in which the non-interacting Green's 
functions are exchanged with interacting ones and vertex cor- 
rections are neglected^. 



Rea xx (u) = ^-^2vo, x (k) 



x / dy np to- nF ( y + u) A{k,y)A{k,v + U ), (6) 



where VQ, x (k) is the bare band velocity in the x direction at 
wave vector k, n F (y) is the Fermi-Dirac distribution function 
and A(k, cj) is the spectral function 



A{k,uj) = ImG K (fc,w) 

TT 



(7) 



where G R (k, uj) is the retarded Green's function. 

Our interest is mostly in low T and low uj properties of 
the conductivity for which the parameter space close to the 
Fermi surface is the most relevant (mainly due to the factor 
(n F (y) — n F (y + u>))/ui). In this parameter space we can 
linearize the bare-band dispersion 

eo(fc) = eo(fcr,0) + v 0iFi r((t))(k F ((f)) - k r ), (8) 

where k F (4>) is a Fermi momentum at angle <fi, which is the 
angle between the (tt, 7r)-(0, tt) and (tt, Tr)-k directions (Fig. 
ID- v o,F,r{4>) is the derivative of the bare band dispersion in 
the k r direction [i.e., the radial from (tt, tt), see Fig. [5). For 
a circular Fermi surface fo,F,r(</>) just corresponds to a Fermi 
velocity. 

By performing the integral over k r the optical conductivity 
can be approximated for a quasi 2D system with 



Re a xx (uj) 



4TT 2 d 



# / dy 



n F (y) -n F (y + uj) 



x —, — Im 



(9) 



d is the distance between CuO layers (d — 11.6 A for 
T12201 35 ), Vq, f (4>) is a Fermi velocity while T, R (c/>, uj) and 
T, A ((j), u) are the retarded and advanced self-energies, respec- 
tively. We assume that they are only ^-dependent in fc— space 
(anisotropic over the Fermi surface) in addition to our pro- 
posed uj and T dependencies. In deriving Eq. (O the integral 
over k r was extended to [— oo,oo] which is a good approxima- 
tion for low T and uj due to the strongly peaked spectral func- 
tion close to the FS. This means any effects of van Hove sin- 
gularities or band edges are neglected. Eq. (0 can be viewed 
as a generalization of Eq. (12) in Ref. |2j to the case of a </> 
dependent self-energy. 

The imaginary part of the optical conductivity can be ob- 
tained by the Kramers-Kronig transformation, 

lma xx (uj) = V dui , (10) 

tt uj'-uj 

or by generalizing Eq. (0 to the complex conductivity 



0~xx(u) 



ie 

A7T 2 d 



d<t> \ dy 



n F (y) -n F (y + uj) 



M^K f(</>) 



v . Fr (4>) uj + -E A (4>,y)-XR(<l>,y + ujy 



(11) 



The plasma frequency uj p is determined by the high frequency 
behaviour (uj 3> band width), 



uj 2 = — lim^oo wlm ct xx (uj) 



(12) 



and in our case this quantity is given by the following integral 
over the Fermi surface 



Air 2 d€r 



^ k F (4>)v 0tF ((j)) 2 

VQ, F , r {4>) 



(13) 



The above expression is equivalent to the band theory expres- 
sion (Ref. Hi (Eq. 13.36)), 



° 2p = vJ ^ nF ^ k ^ 



dhpjk) 
dkl ' 



(14) 



This equivalence can be shown by integrating by parts, con- 
fining the integral to the Fermi surface due to derivative of a 
Fermi function and then using the symmetry a xx = a yy . Here 
eo is the static dielectric constant. 

Using these expressions with our bare band dispersion (see 
Section HUi we obtain uj p ~ 23000 cm" -\ while in Ref. ^ 
they experimentally obtain uj p ~ 15100 cm -1 by integrating 
Rea(w) up to ~ 8000 cm -1 . We believe that this is not a high 
enough frequency to fully exhaust the sum rule. 



A. DC conductivity 

In the limit of uj — > further simplifications can be made, 
1 1 



w + sa^i/J-E^y + w) 2i£"(& y) ' 



(15) 
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where S"(0, y) stands for imaginary part of the retarded self- 
energy. Furthermore, 



riF(y) - riF(y + w) 



dn F (y) 
dy 



(16) 



The DC conductivity can then be written as 



Rea xx 



An 2 d 



d4> 



Vo,F,r(<P) 



dy y -2S"(0,y) 



(17) 



For the bare band dispersion appropriate to T12201 the pre- 



factor 



— - — — — - — turns out to be relatively constant with <fi 
(variation <20%). In comparison the anisotropy of the self- 
energy (1/S" can vary by a factor of more than two) and so 
the pre-factor can therefore be taken out of the integral, re- 
placed with its average value and expressed in terms of uj p 
[Eq. (1131) 1, With this approximation we can rewrite the ex- 
pression for the frequency dependent conductivity in Eq. (flTT i 
in a similar form to Equation (12) in Ref. I22I 



2tt 



dy 



nF{y) -n F (y + u) 



U! 



x dcj> 



(18) 



Using the cos 2 (2(f)) dependence of our model self-energy, one 
can perform the integral over <\> in the Eq. ( TT~8b for to = 0. 
At this point only the integral over the frequency y remains, 
which can be to the lowest order at low T calculated with the 
use of 



dn F (y) 
dy 



(19) 



This is a good approximation, if the self-energy (or 1/S") is 
a fairly constant function of ui for |w| < T, However, further 
improvements can be made by expanding the self-energy part 
to y 2 term and then numerically approximating the y integral 
by Pade approximation, which gives errors less than 10 -6 . 

The resulting expression allows us to perform fits of the 
measured resistivity (p xx = 1/Ke<r xx ) using the three main 
parameters of our model: the strength of impurity scattering 
1/(2to), the strength of AMFL part of self-energy A [where 
\(4>) = A cos 2 (2 </))], and the strength of FL part of self-energy 
s/w* 2 . 

The resulting fits for various T12201 samples with different 
T c s are shown in Fig. [2] Fits to the optimal doping data 
are not performed since they yield unphysical values of the 
parameters (e.g., values of 1/tq ~ 0). This is due to the strong 
increase of the resistivity at optimal doping and is probably 
related to the opening of the pseudogap or some other new 
physics, which is beyond the scope of our model self-energy. 

The resulting fit parameters together with the ones extracted 
from ADMR are shown in Fig. [3] All parameters are con- 
sistent with the ones extracted from ADMR 26 . The zero- 
temperature scattering rate l/(2ro) seems to show an addi- 
tional decreasing trend with increasing T c , which might be 
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Figure 2. Comparison of the measured temperature dependence of 
the intra-layer resistivity p xx to fits for the self-energy model. Fits 
are performed in the T range from above T c to 200 K and agree 
nicely with the data. At higher T > 200 K our model predicts a 
stronger increase in p xx , which could be improved by smoother high 
frequency cutoffs in the self-energy. This would introduce new pa- 
rameters and is beyond the scope of this paper. T c values and refer- 
ences for the data denoted with A and B are given in TableU 



attributed to the loss of interstitial oxygen causing impurity 
scattering. The anisotropic marginal Fermi liquid parameter 
A increases with T c as expected, although it suggests a super- 
linear increase for T c close to the optimal doping. The pa- 
rameter s/uifl is slightly larger than extracted from ADMR 
but still fairly constant with doping. Similar results were also 
obtained from the conductivities of overdoped LSCO 40 . 

Fitting parameters for T c > 70 K become unphysical (too 
small 1/tq) and might be a sign of a new physics out of the 
scope of our simple model self-energy. 

We found that the resulting fit parameter values do not 
change significantly if only the zero frequency self-energy is 
taken into account, as occurs with the delta function approxi- 
mation for the Fermi-function term, Eq. $1% . 

For higher temperatures the measured resistivity shows a 
linear in T dependence over a broader temperature region than 
our model (Fig. [2]). As discussed further below, a smoother 
saturation of the self-energy at high T and high u may im- 
prove the comparison in this regime. 

Saturation of the self-energy may originate in the Mott- 
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Table I. Experimental data sets for the temperature dependence of 
the DC intra-layer resistivity which are fit to our self-energy model. 
(Compare Figures [2] and |3). The corresponding T c 's and references 
are listed. 



P 0.27 0.26 0.25 0.24 0.23 0.22 
0.01 



Data identifier Data T c [K] 



Reference 



A 0, 30, 43, 57, 83 A.W. Tyler et al., Ref. 15 

B 0, 7, 10, 48, 63, 76 T. Manako et al., Ref. .17 
C 15 A.P. Mackenzie et al., Ref. IS 

D 30, 80 A.W. Tyler et al., Ref. 16 



Ioffe-Regel (MIR) limit at which the mean free path I = 
vf,o/(— 2S^ ax ) becomes comparable to the lattice constant 
and electrons become incoherent. Estimate of — S^ ax from 
the MIR limit and our LDA estimate of vf,o ~ la [eV] gives 
— S^ ax ^0.5 eV. This is in good agreement with our maximal 
value of the FL part of self-energy (the main contribution at 
high T) which is ~ 0.5 eV. The MIR limit was already suc- 
cessfully applied to the scattering rate saturation of the opti- 
mal and overdoped cuprates^. It is important to mention, that 
in the underdoped regime, the resistivity saturates at a much 
larger value than expected from the MIR limit, which may be 
due to the smaller carrier concentration^. 



B. Optical Conductivity 

Experiments do not directly measure the frequency depen- 
dent conductivity but rather the reflectivity or absorption of 
a thin film or single crystal. The real and imaginary parts 
of the conductivity are then extracted from a Kramers-Kronig 
analysis^. This is only stable and reliable if there is experi- 
mental data out to sufficiently high frequencies. Furthermore, 
to aid the physical interpretation of the results experimental- 
ists often plot the frequency dependent scattering rate and ef- 
fective mass that is deduced from an extended Drude model 22 . 
However, this also requires a knowledge of the plasma fre- 
quency uj p (compare Eq. (Y2\ ). As mentioned earlier, the bare 
band dispersion from LDA predicts uj v ~ 23000 cm -1 a value 
which is larger than extracted from experiments ( 1 5 1 00 cm - 1 
in Ref. Q. 

To simplify the analysis and avoid the introduction of new 
parameters we compare the results for the model self-energy 
directly to the measured reflectivity. The reflectivity [R(ui)] 
or absorption = 1 — R(u>)] may be written in terms of 

the optical conductivity^ 



R(u) = 

(20) 

in the limit lma xx (u>) equj, which is valid in the fre- 
quency region of the data. Here |o" xa .(ci;)| = ([Rea xx (u))] 2 + 
[Ima xx (u)} 2 y/ 2 . 
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Figure 3. Doping (p) or transition temperature (T c ) dependence of 
the model parameters extracted from fitting the temperature depen- 
dence of the resistivity for a range of T12201 samples. The fitted 
model parameters are consistent with the values and doping depen- 
dencies extracted from ADMR 26 . Strength of the impurity scattering 
1/(2to) shows additional decreasing trend with decreasing doping, 
which might be due to a smaller amount of doped interstitial oxy- 
gen. The AMFL strength A shows a strong increase with decreasing 
doping in good agreement with the results from ADMR^ 6 -. Good 
agreement with the self-energy model is also found for the doping 
independent strength of the FL like scattering s/uj^I. For the doping 
dependence of T c we use the phenomenological relation (see section 
[TTb . T c values and references for the data denoted with A, B, C and D 
are given in Table [TJ 



Comparison of our results, obtained with Eqs. ( f20b and ( [T8b 
and model self-energy parameters extracted from ADMR 26 , 
with the measured absorption is shown in Fig. [4] Agreement 
at low frequencies {oj < 1000 cm -1 ) is quite satisfactory. We 
consider this is quite impressive given that no additional fitting 
parameters beyond those extracted from ADMR 26 have been 
introduced. 

At higher frequencies our model self-energy predicts an ab- 
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Figure 4. Comparison of the measured optical absorption spec- 
tra A(u) for T12201 at three different dopings with the model self- 
energy. The latter is for parameter values extracted from ADMR 26 . 
A good description of the oj, T and T c dependencies of A(oj) is ob- 
tained for oj < 1000 cm -1 . For higher oj, A(oj) shows a stronger 
increase with oj, which could be improved with a softer high-o; cut- 
off for the self-energy. This is similar to what is found for the high-T 
DC resistivity in Fig. [2] Data for T c = 15 K and 70 K are from Ref. 
24 and data for T c = 23 K are from Ref. |H. 



sorption that is too large compared to experimental data. This 
discrepancy could be fixed by incorporating a smoother high- 
frequency saturation making the self-energy more slowly in- 
creasing with oj and rounding its behaviour at the high fre- 
quency cutoff (w^ MFL ~ Ldp L ~ 1700 cm -1 , see Fig. |4|. 
One way of smoothing the high T and to behavior could be 
in adopting the phenomenological approach of Refs. |4l| and 



|44J, where the saturation of the scattering rate is applied by the 
"parallel-resistor" formula, which means the imaginary part 
of the self-energy (Q]) is replaced according to 



1 



1 



1 



s "ff( 



i,T,w) 



yv/ 
^max 



(21) 



Here S 1 'j eal (c6, T, oj) is the self-energy without high-frequency 
cutoffs and S^ ax is the maximal or saturated value of the 
imaginary part of self-energy, and can be treated as a free pa- 
rameter. In the MIR picture described above this parameter 
is estimated to have a value ~ Hvp /a, where a is the lattice 
constant. 

Using two different model self-energies Norman and 
Chubukov 45 performed a detailed analysis of the fre- 
quency dependent conductivity for optimally doped 
Bi2Sr2Ca .92Yo.o8Cu208 +( 5. They deduced a flattening 
of the frequency dependence of the scattering rate near a 
cutoff energy of order 0.3 eV. The high-frequency cutoff may 
also be observed in ARPES as a kink or "waterfall" in the QP 
dispersion due to a noticeable change in dTl jdoj, particularly 
if it obtains a value > 1 (for example see Refs.Hand0). 

The cutoffs give some insight into the underlying physics 
since they tell us the energy scales of the excitations (e.g., spin 
fluctuations, particle-hole excitations) which the electrons are 
scattering offi 4 ^. On the other hand, cutoffs may also reflect 
the limiting scattering rate (e.g., given by the sum rule Eq. 
(0 or the entry into the MIR limit 41 ' 42 ) or entrance into the 
incoherent regime. 



IV. HALL EFFECT 



A. Hall coefficient 



The Hall coefficient in the weak field limit is given by 



Rh 



& xy 



B z [Rea xx } 2, 



(22) 



where B z is the magnetic field in the z or c direction, a xy is 
the Hall conductivity proportional to B z , and Re a xx is the 
in-plane DC conductivity (see Sec. ITTT A) . 

A diagrammatic calculation of the Hall conductivity is 
given in Ref. |46l leading to 



7« - 



-ieB z 



E 



(23) 

where [Ad^B] — Adk^B — (dk^A)B, J M is a current ver- 
tex which we approximate with — evo jfl by neglecting the ver- 
tex corrections, and G R ^ is the retarded (advanced) Green's 
functions. 

The expression in Eq. d23l for the Hall conductivity can be 
further simplified with the following approximations. First, 
we neglect the term <9fc x £'(fc,u;), which arises from differen- 
tiation of the Green's function and is present also as a first 
trivial correction to the vertex, which we also neglect. We 
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find that calculations with this correction do not significantly 
change the results, because our lu) is odd-in-w. Then we 
linearize the dispersion in the k r direction around the Fermi 
surface, Eq. dHJ, and approximate the integral over k r , as was 
similarly done for the DC conductivity (Sec. IIII Ab . Using the 



symmetry a 



(1) - _„W 



Ong in Ref. 47 leads to 



a X y and manipulations similar to those of 



e 3 B z 
4n 2 d 



duj( 



d(j>[-Vo fF {4>) x d^vo^ (</>)} 
dnp(uj) 1 



dtu ; (-2£"(0,w)) 2 



(24) 



where we have also used that our self-energy depends only 
on <fi in momentum space. A more detailed derivation can be 
found in Appendix IB! 



B. Comparison with the Boltzmann equation 

Ong has given an elegant geometrical interpretation of the 
Hall conductivity a xy for a two-dimensional Fermi liquid 47 . 
It is proportional to the area swept out by the scattering length 
or mean-free path 1(0) = vf(<^)t(0) as one traverses the 
Fermi surface. This illustrates how the Hall effect is sensi- 
tive to anisotropy in the Fermi surface via the Fermi velocity 
vp(</j) and and to anisotropy in the scattering time r(4>). 

Eq. (l24l is consistent with the expression derived from the 
Boltzmann equation 41 and with Ong's geometric expression 47 . 
If the frequency dependence of the self-energy close to lu = 
is neglected, then 



dn F (u) 

M — h— ) 



1 



1 



dcu '(-2E"(0,a;)) 2 (-2S"(0,O)) 2 ' 

(25) 

To make the comparison with the Boltzmann equation and 
relaxation time approximation more explicit, we start with the 
Boltzmann equation result for the Hall conductivity 4 -^, 



2n 2 d 



J d 2 k{- 



^f^)^vxB.V(^). (26) 



Here Ek is the quasi-particle (QP) dispersion, v = (v x ,v y ) 
is the QP velocity, and V is the scattering rate, which are all 
fc-dependent. The integral goes over the first Brillouin zone in 
two dimensions. Symmetrizing the expression with the use of 



a^x and applying 



|vxB.V(^)=vxB.(iv(^)-^V(o),. ,27> 



leads to the following expression, 

= 3 r dn F (E k 



a 



(i) 



4n 2 d 



d A k(- 



dE k 



W x B 



•(«*V(t>„) -v y V(v x ))^. (28) 



Furthermore, if the integral over the 2D Brillouin zone is de- 
composed into the integrals over <f> and k r and in addition the 



QP dispersion is linearized close to the Fermi surface with 
Ek ~ VF ir (<f>)(kF(<f>) — k r ) and Y ~ Y(<j)), then the integral 
over k r may be performed (neglecting band edge effects) and 
we are left only with the integral over <\>. For the magnetic field 
in the z direction one can then rewrite a2y in a similar form 
as in Eq. d24T i if the integral over lu in Eq. (l24b is replaced 
with t(a\z (similar to Eq. d25b). We should note here that 
the expression derived from the Boltzmann equation includes 
only renormalized quasi-particle entities, while Eq. (l24l > in- 



cludes only non-renormalized quantities. This is not a prob- 
lem since the renormalization cancels by taking vp = Zvo,f 
and r = — 2ZS". However, this might not be the case, if the 
shape of the non-interacting Fermi surface is changed due to 
the renormalization. This does not happen for our model self- 
energy, since its real part is always zero at to = due to the 
imaginary part being an even function of frequency. 

The relationship to Ong's geometric interpretation is more 
straight forward. If the integral over u> in Eq. d24l is approxi- 
mated as in Eq. d25l ), one can write 



47T 2 d 



d#(0) X fyl(^)], 



(29) 



where 1(0) = Vo } p((j))/[— 2E"(</>, 0)] is the mean free path 
used in the Ong's geometrical interpretation of the Hall con- 
ductivity. From this expression it is nicely seen that the renor- 
malization Z cancels and that the Hall conductivity is deter- 
mined by the mean free path on the Fermi surface. 



C. Comparison with experiment 

The zero temperature (T = 0) value of the Hall coeffi- 
cient Rh for a circular Fermi surface corresponds to l/(en e ) 
with n e being the density of electrons. Deviations from this 
value depend on the shape of the Fermi surface. If for our 
model tight-binding band structure we assume a rigid band 
shift from the highly overdoped to optimally doped regime, 
then the T — value of Rh is expected to change by less 
than 10%. Temperature broadening affects Kea xx and a^y 
only at higher T and is within our model estimated to result in 
a relative change of - 10% at T ~ 200 K. The T broadening 
effect is reduced in Rh and is estimated to be < 0.5 x 10~ 10 
m 3 /C. In contrast to the above relatively small changes with 
temperature and doping, experiment shows that Rh(T) can 
vary by as much as 100%, as T c varies from to 50 K in the 
overdoped regime (compare Figure |7). 

As can be seen from Eq. d24"l > for a x X y and Eq. ( TPTI i for 
Kea xx the temperature dependence of the Hall coefficient 
comes from the T-dependence of the anisotropy of the scat- 
tering rate. This becomes more apparent if we rewrite the Hall 
coefficient in the following form 



R 



(30) 



where we have neglected the T-broadening effects, fnii 1 ) 
is the (f> dependent coefficient (corresponding to the Hall 
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conductivity), which needs to be integrated over <fi and de- 
pends only on bare-band dispersion. /dc(0) i s similar to 
fn (</>), but for the DC conductivity (see Appendix [A}. The 
only T-dependent quantity in the above equation is the self- 
energy and its T-dependent anisotropy is responsible for In- 
dependent Rh- This is in agreement with results in Ref. HH 
However, the absolute change of Rh with temperature de- 
pends strongly on the shape of the Fermi surface. This is 
demonstrated in Figs. |5]and|6] 



(71, 




1 X W 



Figure 5. (Left) Three different approximations for the shape of the 
Fermi surface (FS). The FS deduced from ADMR measurements 12 , 
ARPES measurements 8 , and LDA calculations^ 7 , are denoted with 
ADMR, ARPES, and LDA, respectively. (Right) Curves of the 
mean-free path 1(0) as one goes around the FS for the three dif- 
ferent FSs. According to Ong's geometric interpretation 47 the encir- 
cled area is proportional to a^y [Eq. I129H. Although the shapes 
of the FSs do not change much between different approximations, 
the mean-free paths and the encircled areas in 1 space change sub- 
stantially. The main difference comes from the curvature of the 
FS (90Vo,f(0)) close to <f> ~ tt/8. The absolute value of Rh is 
therefore very sensitive to the shape of the FS. The mean-free path 
1(0) was calculated with our model self-energy for T = 100 K, and 
T C = 30K. 

The overall doping (or T c ) and T dependence of the mea- 
sured and calculated Rh are shown in Fig. [7] and [8] respec- 
tively. The temperature dependence of Rh suggests, that the 
scattering anisotropy strongly (linearly) increases at low T (in 
our model due to the AMFL part of self-energy), reaches its 
maximum at ~ 110 K and then the scattering slowly becomes 
more isotropic again as the FL part of the self-energy model 
begins to dominate. 

The fact that for T c = the experimental Rh shows a 
small T-dependence (see Fig. Q represents a problem for our 
model, since the model has no anisotropy for T c = 0. How- 
ever, there was no ADMR data for T c = and so it is possible 
that the anisotropy actually does not go to as T c 0, or per- 
haps that our assumption that the T 2 term is strictly isotropic 
needs to be relaxed. 

Comparison of our results (Fig. [8) with the measured Rh 
(Fig. 13 shows qualitative, and to some extent also quantita- 
tive, agreement. However, our Rh does reach a maximum for 
T ~ 80 K, while the maximum appears at higher T (~ 110 K) 
in experiment (Fig. |7). In order to get a better comparison the 
FL part of our self-energy model should be reduced (smaller 
s/ujfl). Also inclusion of a smoother high frequency cutoff 
could move the maximum in our Rh (T) to higher T. 
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Figure 6. Calculated temperature dependence of the Hall coefficient 
(Rh) for three slightly different Fermi surfaces (see Fig. [3}. The 
absolute change or maximal value of Rh depends strongly on the 
shape of the Fermi surface. All curves are calculated with the same 
scattering rate, using our self-energy model for T c = 30 K. Temper- 
ature broadening effects due to the Fermi-Dirac distribution are not 
taken into account since they are small (< 10%). The slightly dif- 
ferent values of Rh at T = 0, where the scattering is dominated by 
impurities and is therefore isotropic, also comes from small changes 
in the Fermi surface shape. 




T[K] 

Figure 7. Temperature dependence of the measured Hall coeffi- 
cient Rh for several different T c s. These all correspond to dop- 
ings for which there is a large hole Fermi surface and show a non- 
monotonic temperature dependence that increases with increasing 
T c . For T < T c Rh may be strongly suppressed by the supercon- 
ducting transition. Data for T c = K and 8 1 K are for polycrystalline 
samples measured in Ref. 20, data for T c = 10 K and 50 K are from 
Ref .Til T = 15 K data is from Ref. QJ, and T c = 25 K data is from 
Ref.lJ 



In fitting our model to Rh it turns out that one parameter is 
free (one of 1/ (2to), A or s /u>£i). This is because the absolute 
value of Rh is unchanged by a re-scaling of the scattering 
time, as can be seen from Eq. d30l ). This is closely related to 
Rh not depending on r in a simple FL picture. 
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Figure 8. Temperature dependence of the Hall coefficient Rh calcu- 
lated with our model self-energy for several T c s and for the ADMR 
Fermi surface 12 . Results should be compared with Fig. [7] The figure 
illustrates how decreasing doping (increasing T c or A or anisotropy 
of the self-energy) leads to a large change in the magnitude of the 
temperature dependence. All results are obtained for the same Fermi 
surface and only the anisotropy of the self-energy is changed. Ar- 
rows shown in the lower right indicate the weak dependence on the 
band structure change with doping, indicating the absolute shift of 
the zero temperature value of Rh (T c = 0) for two different dopings 
(p = 0.15 and p — 0.30) as given by a rigid band shift. 



D. Hall angle 



The Hall angle is defined by 

Rea xx 



cot 9 H (T) 



Pxx{T) 

Rh(T)B z 



(31) 



Since our model can describe the temperature dependence of 
the intra-plane resistivity and Hall coefficient, as we showed 
above, it must also describe the Hall angle. Here, we exam- 
ine the temperature dependence of cot Oh in order to point 
out that the observed T 2 dependence of cot Oh (cf. Figure |9]l 
naturally follows from our model self-energy and that there is 
therefore no need to evoke more exotic theories in order to ex- 
plain the qualitatively distinct temperature dependence of p xx 
and cot Oh- 

Experimental data and our results are shown in Fig. [9] and 
provide additional support for our model self-energy. In par- 
ticular, the linear dependence of cot Oh on T 2 supports the 
T 2 dependence of the isotropic part of self-energy or scat- 
tering rate in the nodal direction. That is because cot Oh is 
dominated by the isotropic part (Sp L ), while it suppresses 
the anisotropic part (S^mfl) of the self-energy. This point 
was previously emphasized by Carrington et al.,— Ioffe and 
Millis 50 and by Stojkovic and Pines 48 (see also Ref. l5lt) . To 
show this more explicitly, we use a similar expression to the 
one in Eq. (|30l ), approximate /#(</>) an d fDc( ( t ) ) with a con- 
stant, and perform the integrals over <f>. This leads to 



cot Oh oc 



-£fl(0)[1 



££l(0) 



£W(o,o) J 



(32) 



It turns out that the temperature dependence of cot Oh is dom- 
inated by the first factor, because the second factor is weakly 
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Figure 9. Temperature dependence of the Hall angle cot Oh- (Top) 
Experimental data show that cot 6h has a linear dependence on T 2 
and a weak doping dependence. A similar dependence is found for 
our model (Bottom), where the linear dependence of cot 8h on T 2 
comes predominantly from the Fermi liquid like scattering (Sp L ) that 
gives the scattering rate on the nodal part of the Fermi surface (which 
also gives the dominant contribution to the Hall conductivity, see Fig. 
[5}. (Top) The Inset shows a small down-turned deviation from the T 2 
dependence of cot 6h at low T (for T c = 15 K). This is also obtained 
within our model (bottom inset), although not as pronounced as in the 
experimental data. Experimental data for T c — 10 K and 50 K are 
from Ref. 
is from Ref. 



T c = 15 K data is from Ref. 0J and T c = 25 K data 
2U Our model results are calculated with the ADMR 



Fermi surface. 



temperature dependent. For more details see Appendix [C] 
Hence, the Hall angle is dominated by isotropic scattering or 
by the region on the Fermi surface with the weakest scattering 
or the longest mean-free-path, while the effect of anisotropic 
scattering is suppressed. Further suppression of anisotropic 
part comes from the anisotropy of which is larger in 

the nodal and smaller in the antinodal direction. 

Although the effect of S^mfl on cot Oh is small (note T 2 
dependence in Fig. [9]), it still changes the pure T 2 dependence 
of cot Oh to T n with n < 2. Values of n < 2 were actually 
observed in YBCO (Ref. 52) and Bi2201 (Ref. HI and Hi) 
where n changes from ~ 1.8 in the optimal or underdoped 
regime to n ~ 1.6 in the overdoped regime. Our model pre- 
dicts n = 2 in the highly overdoped regime where the AMFL 
part of the self-energy is zero, but could predict n < 2, if the 
smoother high-frequency cutoff were introduced. This would 
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make the T 2 dependence of FL like self-energy more linear in 
T for higher T, observed experimentally in Bi2201 (Ref. l4lt) . 
However, with decreasing doping and consequently increas- 
ing anisotropy our model would predict a further decrease of 
n, which is the opposite trend to that observed experimen- 
tally Bi220 1 51 ' 53 . In contrast, a different model with strong 
anisotropic impurity scattering, an anisotropic term oc T 2 , and 
a smooth high-frequency saturation yields an increase of n 
with increasing anisotropy 41 . For T12201 no change of n with 
doping was observed so far, which might be due to a more 
square-like Fermi surface and therefore the decreased effect 
of anisotropy on cot Oh- 

Our anisotropic self-energy model is therefore capable of 
simultaneously describing the linear in T part of the DC con- 
ductivity and T 2 dependence of cot Oh over a wide doping 
range from optimal to the heavily overdoped region. This 
shows, that there is no need to introduce more exotic theories 
with two types of quasi-particles (e.g., spinons and holons) 
with different scattering rate a 54,55 , to capture the qualitatively 
different temperature dependence of p xx and cot Oh- 



V. COMPARISON WITH MICROSCOPIC MODELS 

It is a challenge for microscopic theory to quantitatively 
describe the observed temperature and doping dependence of 
transport properties or equivalently the proper T, p and <f> de- 
pendence of the self-energy in the overdoped cuprates. In this 
section we compare our model self-energy to results from sev- 
eral microscopic theories in order to evaluate their potential 
for a successful description of the various experimental data. 

A weak coupling treatment of the Hubbard model can pro- 
duce an anisotropic scattering rate of similar frequency and 
angular dependence to our model. The anisotropic MFL 
component arises from a nesting of the Fermi surface in 
the anti-nodal regions 56 or from proximity to a van Hove 
singularity 56 ' 57 . However, for the latter case the resulting scat- 
tering rate would have the opposite doping dependence and 
would appear only at a higher temperature than that experi- 
mentally observed for T12201, since the van Hove singularity 
would reach the Fermi surface for dopings larger than in the 
highly overdoped regime. Hence, the anisotropic MFL term 
can only arise from the nested parts of the Fermi surface which 
produce a particle-hole susceptibility similar to that found in 
one dimension. Hence, the scattering is essentially arising 
from particle-hole excitations with a high-frequency cutoff of 
the order of the band width. 

A functional renormalization group treatment of the Hub- 
bard model 58 (for a review see Ref. l59h shows a T, p and 
(f> dependence of the scattering rate in qualitative agreement 
with ADMR 1112 and our self-energy model. However, it pre- 
dicts an order of magnitude smaller anisotropic scattering rate 
than observed in experiment, while it gives the correct order 
of magnitude for the isotropic scattering (oc T 2 ) in agreement 
with our self-energy model (see supplemental material of Ref. 

HI). 

The Hidden Fermi liquid (HFL) theory by Casey and 
Anderson 60,61 uses a Gutzwiller projection of the Fermi liq- 



uid wave function. However the scattering rate predicted by 
HFL has a linear T dependence only for temperatures above 
T ~ 400 K, in strong contrast to the ADMR measurements 11 , 
where the T linear term is observed even for T < 60 K 26 . 
Furthermore, within the HFL theory the anisotropic scatter- 
ing emerges solely as a consequence of anisotropy of the 
Fermi momentum and of the Fermi velocity on the Fermi 
surface 61 - 62 . LDA calculations 37 show a weaker anisotropy 
and with the opposite doping dependence than that needed in 
HFL to capture the experimentally observed scattering rates 26 . 

Cluster dynamical mean field theory (CDMFT) can also 
calculate scattering rates at different parts of the Fermi sur- 
face. Results presented in Ref. 63 and obtained with a Hub- 
bard model with t' jt = —0.15 and U — It reveal quali- 
tatively similar behaviour to ADMR and to our model self- 
energy. For higher dopings CDMFT gives an isotropic scat- 
tering rate, which becomes more anisotropic (stronger scatter- 
ing in antinodal direction) and stronger with decreasing dop- 
ing. However, due to limitations of the quantum Monte Carlo 
method CDMFT is currently limited to T > 0.05i ~ 200 
K, which is above the most interesting experimental regime. 
Quantitative comparison with our self-energy model shows, 
that CDMFT 63 predicts at T = 200 K a smaller isotropic 
part, by a factor ~ 2.5. Comparison of the anisotropic part 
is complicated due to patch averaging in DMFT. However, the 
CDMFT self-energy 63 has the same order of magnitude as our 
model self-energy, at least at T ~ 200 K. Detailed quantita- 
tive comparison with the CDMFT results is given in Appendix 

is 

Treatment of the t-J model with the finite-temperature 
Lanczos method (FTLM) 6 ^ yield results in good agreement 
with several experimental data, including the optical conduc- 
tivity and high T resistivity. However, the temperature range 
of reliable results (due to finite size effects) obtained with 
FTLM is too high to address the low T transport properties 
and in particular the anisotropy in the scattering rate observed 
in ADMR. 

A large-iV expansion treatment of the t-J model 65 , found a 
scattering rate with a similar temperature and angular depen- 
dence as our model self-energy. However, as optimal doping 
is approached it also exhibits divergence of the anisotropic 
scattering rate at low temperature, due to a rf-density wave in- 
stability near optimal doping. This is qualitatively different 
from our model self-energy. 

Ioffe and Millis 50 considered how superconducting fluctua- 
tions could produce an anisotropic scattering rate. They sug- 
gested that in the overdoped region the rate should scale with 
T 2 , but it should be kept in mind this depends on what as- 
sumptions one makes about the temperature dependence and 
magnitude of the superconducting correlation length. 

Metzner and colleagues have been investigating d-density 
wave fluctuations near a quantum critical point associated 
with a Pomeranchuk instability. 66 Their starting point was an 
effective Hamiltonian which has a d-wave form factor built 
into it. But this was motivated by earlier work 67 on the Hub- 
bard model which found from renormalisation group flows 
that strong forward scattering developed led to a Pomeranchuk 
instability. Although, this work reported an anisotropic scat- 



11 



tering rate that is linear in temperature it turns out that due to 
vertex corrections the transport scattering time scales as T 4/3 
and the resistivity scales as T 5/3 (Ref. roa) . 

In spite of all these theoretical studies the question remains 
whether there is a simple explanation for the scattering in 
terms of a single mechanism: e.g., antiferromagnetic, super- 
conducting, or d-density wave fluctuations. Furthermore, is 
there a smoking gun experiment which could distinguish be- 
tween these different contributions? For example, they should 
have a different dependence on the magnitude of an external 
magnetic field. We also note that a magnetic field couples dif- 
ferently to spin and orbital degrees of freedom, and the former 
contribution is dominant for fields parallel to the layers. 



VI. CONCLUSIONS 

In conclusion, we have shown that our model self-energy 
can describe a wide range of experimental data on overdoped 
cuprates. In earlier work we showed it could describe scat- 
tering rates deduced from ADMR, the quasi-particle disper- 
sion seen in ARPES, and effective masses deduced from spe- 
cific heat and quantum magnetic oscillations 26 . Here, we have 
shown that neglecting vertex corrections the model can also 
describe experimental data on electrical transport properties, 
including DC conductivity, optical conductivity, Hall coeffi- 
cient, and Hall angle. 

The small quantitative discrepancies between the model 
and measured data at high frequencies (>1000 cm -1 ) or 
higher T (>300 K) could be reduced with application of a 
smoother high-frequency cutoff for the self-energy, e.g. with 
the "parallel resistor" formula 41 ^. 

The successful description of the experimental data by our 
analysis shows that inclusion of vertex corrections is not nec- 
essary at this level of approximation. However, for the Hub- 
bard model on a square lattice it is claimed 69 70 that vertex cor- 
rections are important in the optimal and underdoped regimes. 

Our results on the DC resistivity show that in the overdoped 
regime the isotropic scattering weakly depends on doping (or 
T c ), while the anisotropic scattering increases super-linearly 
with increasing T c of decreasing doping. Similar findings 
were obtained for LSCO in Ref. |40j This highlights the fact 
that the doping dependence of the DC resistivity in cuprates is 
generic and not so dependent on material properties or Fermi- 
surface shape. 

Such generic behaviour is not seen in the Hall effect, where 
for overdoped LSCO the Hall coefficient monotonically de- 
creases with increasing temperature and increasing doping 71 , 
with a sign change for a doping p ~ 0.3. This may be due to 
the proximity of the Fermi energy to the van Hove singularity 
in LSCO. 

We have also shown that the main temperature dependence 
of the Hall coefficient Rh comes from the temperature depen- 
dence of the self-energy anisotropy. Our model was contrasted 
with the Marginal Fermi liquid (MFL) model of Abrahams 
and Varma 28 , which consists of an anisotropic impurity scat- 
tering term and an isotropic marginal Fermi liquid term. This 
model was used to describe the T dependence of the Hall an- 



gle at optimal doping. However, their model cannot describe 
the pronounced non-monotonic T dependence of Rh found 
in overdoped T12201. It may be worth noting, that overdoped 
LSCO, in contrast to T12201, shows a monotonic T depen- 
dence of Rh and so may be adequately described by the MFL 
model 2 ^. 

On the other hand, the observed T n dependence with n < 2 
of cotOn is generic in the cuprates and has in combina- 
tion with the T-linear resistivity stimulated the proposal of 
more involved theories. For example Anderson 54 suggested 
two types of quasi-particles with different scattering rates. It 
was suggested that, different scattering mechanism may be 
connected to the charge conjugation properties of different 
currents 55 . However, our analysis shows, that there is no need 
to evoke such theories, since our anisotropic self-energy gives 
consistent quantitative description of both p xx and cot 6h- In 
addition, we have shown that it also quantitatively describes 
the frequency dependent conductivity, remarkably with no ad- 
ditional fitting parameters and just using the parameters orig- 
inally extracted from ADMR2 6 .. 

Future work would could and should consider calculation 
of thermoelectric transport properties such as the Seebeck co- 
efficient and Nernst signal using the same model self-energy. 
In a quasi-particle picture both of these transport coefficients 
contain contributions from the energy dependence of the scat- 
tering time 72 ' 73 and so may be sensitive to a marginal Fermi 
liquid contribution to the self-energy. 

The relevance of the model self-energy to electron doped 
cuprates 74 should also be investigated. Recently it was 
observed 75 that in the overdoped region of the phase diagram 
the resistivity had a linear-in-temperature term which was pro- 
portional to the superconducting T c , as in the hole doped 
cuprates considered here. 

The broader significance of this work is that it shows that 
the metallic state in the overdoped regime is not a simple 
Fermi liquid and exhibits some physics which is similar to that 
found at optimal doping [marginal Fermi liquid behaviour] 
and underdoping [anisotropic Fermi surface properties with 
cold spots in the nodal directions]. A significant challenge is 
to find a general phenomenological form of the self-energy 
that with decreasing doping smoothly crosses over to a form 
that describes the pseudogap state, such as the form proposed 
by Yang, Zhang, and Rice 5 ^ 6 .. 

ACKNOWLEDGMENTS 

This work was supported by an Australian Research Coun- 
cil Discovery Project (DP1094395) and the EPSRC (UK). We 
thank K. Haule, J. Merino, P. Prelovsek, B.J. Powell, and J. 
Schmalian for helpful discussions. NEH would also like to 
acknowledge a Royal Society Wolfson Research Merit Award. 

Appendix A: Functions fn (<f>) and fDc(<t>) 

Here we give explicit forms for the functions Jh{4 > ) an d 
/uc(^) that appear in Eq. ( f30l >. The function can be 
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readily obtained from Eq. d24i i for a X y , 



M<f>) 



A-K 2 d 



(Al) 



On the other hand, fDc(4>) can be obtained from Eq. (l% for 

Re a xx , 



fnc(<t>) 



e 2 M^K.fM 
Air 2 d VQ,F, r {4>) 



(A2) 



Appendix B: Details on derivation of Hall conductivity 

Here we give more details of the derivation of the expres- 
sion for the Hall conductivity, Eq. (l24l l. by starting with Eq. 
( 1231 , which is taken from Eqs. (2.7) and (3.36) in Ref. |46j J x 
in Eq. (1231 represents the current vertex, which in our approx- 
imation of neglecting vertex corrections equals — ev x . The 
square brackets denote 



[Ad fl B] = Ad ku B-(d k A)B, 



which leads to 



[J X dyJy] = e 2 (V x d ky Vy ~ Vyd ky V x ), 



(Bl) 



(B2) 



[G R d x G A ] = G R d kx G A ~ d kx {G R )G A 
= G R G A {G A {hv x +d k Y A ) 

-G R (hv x +d k X R )). (B3) 

qR(A) re p resen t the retarded (advanced) Green's function, 
which may be written in terms of retarded (advanced) self- 
energy With the use of T, R — Y, A * = E = E' + iS" 
and G R ^ =l/(u-e k - E fl ( A )) we can write 

[GR ~ d * GA] = [(.-e k Jr + (^ (B4) 

x [(hv x + <9 fcx E')(-2iE") + (-2id kx X")(u - e k - £')] • 
The Hall conductivity can now be written as 



'xy 



-ie 3 B z ^ f dui dn F (ui) . 
O / ; / <->_ V 



2ir 



doj 



(B5) 



x [{v x + ^E')(-2iE") + {-2id k Y"){u -e k - E')] . 

Since we neglected vertex corrections we should also neglect 
d kx E' in the above equation, which is the same as neglecting 
the first correction to the vertex. For our even-in-w E" first 
order vertex corrections (v x — > v x + d km E') turn out to be 
negligible. Also the term with (oj — e k — E') may be neglected 
due to the strongly peaked and even-in-w prefactor l/[(w — 
e k - E') 2 + (E") 2 ] 2 . 



r(l) 



e A B 
~~ 2 

xv(k) 



-E 



doj dnp(u)) 
-2E" 



- e k - E') 2 + (E") 



(B6) 



where 



)(k) = v 2 x d ky v y - v x v y d ky v x 



(B7) 



The sum over k may be converted to an integral over the first 
BZ, and the integral over k z can be performed due to the 
quasi-two-dimensional nature of the system. The integral over 
k x and k y may be decomposed into integrals over k r [the ra- 
dial direction from (it, it), see Fig. [5) and its azimuthal angle 
6. We are left with 



xv {2n) 3 d 

x / dk r k r v(k) 



-2E" 



[(oj - e k - E') 2 + (E") 2 ] 2 



.(B8) 



In the next step we linearize the bare band dispersion close to 
the Fermi surface in the k r direction [see Eq. ©I and approx- 
imate 



1 



[(w - e k - E') 2 + (E") 5 



2(-E"(</>, W ))V^W>) 



<5[fc r - k r {<f>,u% (B9) 



with 



k r ((f),uj) = k F (<f)) + 



LJ — S'(0, w) 
^O,F,r(0) 



(BIO) 



We further approximate k r (<p,uj) ~ k F ((f>), which we have 
checked numerically results in an error of less than 2% for the 
relevant band structures. With this approximation the integral 
over k r can be explicitly evaluated. 



e 3 B x 



do 



k F {4>)v(k F {4>),4>) 



(1) _ 

-~d J ^ Vo,F, r (<f>) 

dn F {u) 1 



x / du(- 



duj / (-2E"(</),w)) 5 



(Bll) 



There is one further simplification regarding the "velocity" 
term that can be done. Using the symmetry a x J = —cr yx 
we can write Eq. dB7b 

2 ~^ \\ v tdk y v y + Vyd kx v x - v v v x (d ky (v x ) + d km (v y ))\ 

= x e z) ■ (v v Vv x - v x Vvy), (B12) 

Expressing (v x e z ) = vt, where t is unit vector parallel to 
the Fermi surface, and using 



Vt>, • t = 



Vv x ■ tdk\\ _ d^,v x {4>) 



dkw 



fc/(0)/cos< 



(B13) 



where 9 is the angle between the Fermi surface direction and 
direction (perpendicular to k r ). Analysing in the same way 
the y term brings us to 



1 / y \ r a COS6» 

V -> -«O,F(0) -VO,F X d^Vo.FWz , /, v 

2 w 9) 



(B14) 



Finally, using Uo,F,r(0) = cos0uo.f(0) cancels cos 6* and we 
can write our result as Eq. d24l i. 
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Appendix C: Effect of anisotropy on Hall effect 

Here we demonstrate how the anisotropy in the scattering 
rate influences the Hall effect. In particular, we show with 
a simple example that the T-dependence of the Hall coeffi- 
cient Rh is dominated by T-dependent anisotropy, while, on 
the other hand, the Hall angle cot Oh and its T-dependence 
are dominated by the isotropic scattering. We start with the 
expressions for conductivities, which were used in obtaining 
Eq. ([3(3, 



= / d4>f DC {4>)- 



l 



-2£"O,0)' 



(-2£"(</>,0)) 2 



(CI) 



(C2) 



In this simple approximation we neglect the (f> dependence of 
functions fni't') an d /dc(0) an d exchange them with then- 
average values fn and foe- This is feasible due to the much 
stronger anisotropy in the self-energy than in the /-functions. 
Further on, we use a shorter notation for the two self-energy 
parts, — Sp L (0) = a and — Eamfl(0, 0) = b, which allows us 
to write 



£"O,0) = a + 6cos 2 (2(/0, 



(C3) 



where a and b are T dependent, a includes impurity scattering 
and the FL like part which is oc T 2 , while b is due to the 
AMFL part and is oc T. With this approximation, integrals 
over 4> in Eq. dCU and dC2b can be explicitly performed and 
lead to 



1 



DC 



a(l + |)Va' 



7T - 2 + ^ 

4 ^2(1 + |)3/2 ' 



(C4) 



(C5) 



Expressing the Hall coefficient and Hall angle in this approx- 
imation brings us to the final result of this section, 



Rh = 



4tt/ 2 



DC 



"(I 



1 



l + b/a 



(C6) 



or decrease the contribution from the AMFL part of the self- 
energy (or the antinodal part of the Fermi surface). The ef- 
fect of changing the /-functions by changing the shape of the 
Fermi surface can for example be seen in Fig. [6] Furthermore, 
the effect of the anisotropy factor (1 + 1+ ^y a )~ 1 is to down- 
turn the T 2 dependence of cot Oh and make it more like T n 
with n < 2, which has in fact been observed (see inset in Fig. 

mandRefs.EUandHi) 



+ 0.5 



T c = OK 

T c =20 K 

T c =40 K 

T c =60 K 

T =80 K 
. i 



100 



200 



300 



T[K] 



Figure 

l + b/a 



10. Small temperature dependence of the factor (1 + 
_1 appearing in cot 9h- This factor shows less than 10% 



variation with T from 50 to 300 K. Therefore the main T depen- 
dence of cot 6h comes from the isotropic scattering, which in our 
model self-energy is oc T 2 (in agreement with experiment). 



Appendix D: Comparison of model self -energy with CDMFT 

Here we show a quantitative comparison of our model self- 
energy with Cluster Dynamical Mean-Field Theory (CDMFT) 
calculations on the Hubbard model 63 . Scattering at the Fermi 
surface or £"(w = 0) is the most relevant quantity for expla- 
nation of many transport data, which we analyze in this work. 
Our model = 0) at a doping level p = 0.3 is com- 

pared with CDMFT results in Fig. [TTJ Comparison of the de- 
pendence of the self-energy on the Matsubara frequencies on 
imaginary axis (see Fig. H~2b can be done to avoid analytical 
continuation of CDMFT results. The slope at low frequencies 
(9 Wn Im£(iw n )| u?i ^o) is related to the quasi-particle weight 
and mass renormalisation. 



cot I 



'H 



JH 



1 



l + b/a 



(C7) 



From Eq. iC6i it is evident that the Hall coefficient, and in par- 
ticular its T dependence, are dominated by the T-dependent 
anisotropy b/a. On the other hand, Eq. (IC7l i reveals that 
the Hall angle cot Oh is dominated by the isotropic scatter- 
ing a, while the anisotropy effect is strongly suppressed in 
the factor (1 + 1+ ^y a )~ 1 . The doping and T dependence of 

(1 + 1+ 1 b/ / a ) -1 for our model self-energy are shown in Fig. 
[Tol The effect of anisotropy can be further increased or de- 
creased by <j> dependent Jh or /dcs which can either increase 
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T[K] 



Figure 11. Comparison of the imaginary part of the self-energy 
at zero frequency (-E"(w = 0)) from CDMFT results (Ref. 
Fig. 11) with our model self-energy. CDMFT results were obtained 
for different doping levels p — 0.3, 0.25, and 0.2 and for different 
patches on the Fermi surface ("nodal" denotes nodal patch, while 
"antin." denotes antinodal patch). CDMFT results are only avail- 
able at higher temperatures (T > 200JO due to limitations of the 
quantum Monte Carlo method used. Our model self-energy is most 
reliable at low T. (Experiments suggest it should become more lin- 
ear for T > 200 K (see Fig. [2}). Most reliable comparison with 
CDMFT can therefore be done at T = 200 K. At such T CDMFT 
predicts a weaker isotropic self-energy (compare CDMFT n = 0.70 
and our nodal self-energy). Quantitative comparison for stronger 
anisotropies or lower doping is more difficult due to patch averag- 
ing in CDMFT. However, CDMFT predicts the correct trend with 
doping and order of the magnitude for the self-energy. Our antinodal 
self-energy was calculated with T c — 60 K, and the energy scale of 
CDMFT data was set with the LDA hopping parameter ti = 0.438 
eV. 



Figure 12. Comparison of results for the Matsubara frequency 
dependence of ImE(iw n ) from CDMFT (Ref. 63 Fig. 8) with 
our model self-energy. It is seen that the saturated value of the 
self-energy is similar in both results. The isotopic CDMFT result 
(p = 0.28) predicts a smaller slope at low frequencies (and thus are 
a quasi-particle weight closer to one) than our model. Parameters for 
our model self-energy are the same as in Fig. 1111 and the tempera- 
ture corresponds to that of the CDMFT calculation, T — 0.05£i ~ 
250K. 
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